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Abstract 

Motivated by isotropic fully developed turbulence, we define a theory of symmetric 
matrix valued isotropic Gaussian multiplicative chaos. Our construction extends the 
scalar theory developed by J. P. Kahane in 1985. 

1. Introduction 

In the pioneering work [15], J. -P. Kahane introduced the theory of Gaussian multiplicative 
chaos. Given a metric space and a reference measure, Gaussian multiplicative chaos gives a 
mathematically rigorous definition to random measures defined as limits of measures with 
a lognormal density. The main application of this theory was to define the Kolmogorov- 
Obhukov model of energy dissipation in a turbulent flow (see [17], [19]): in this context, 
the metric space is the euclidean space M'^ equipped with the Lebesgue measure and the 
log density has logarithmic correlations. Since this seminal work, the theory of Gaussian 
multiplicative chaos has found many applications in a broad number of fields among which 
finance ([2], [10]) and 2-d quantum gravity (see [9], [18] for the physics literature and [3], 
[11], [22] for the mathematics literature). 

The main motivation of the Kolmogorov K41 theory ([16]) and it's extensions ([17]) 
is to define a realistic statistical theory of an incompressible, homogeneous, isotropic and 
fully developed turbulent flow (see for example [12, 13]). This ambitious program consists 
in defining a probabilistic model for the velocity field which satisfies the main statistical 
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signatures observed experimentally, such as the mean energy transfer towards the small 
scales and the intermittency (or multifractal) phenomenon ([16]). Ideally, one looks for a 
field as close as possible to an invariant measure of the equations of motion. In [6], the 
authors propose a probabilistic construction of such a velocity field. Their construction, 
which requires a limiting procedure, is mathematically non rigorous and is based on the 
short time dynamics of the Euler flow, as well as further multifractal considerations. One 
of the key step of this construction is the introduction of the exponential of an isotropic 
trace-free matrix whose entries are Gaussian variables with logarithmic correlations. Un- 
fortunately, as far as we know, there is no matrix valued theory of Gaussian multiplicative 
chaos. The purpose of this work is thus to define such a theory for a class of Gaussian 
symmetric and isotropic matrices. 

In the next section, we present the framework and the main results. Section 3 is 
devoted to the proofs of our main results. In the appendix, we gather general formulas 
which are useful in our proofs. 

Notations: we denote by M{W^) the set of measures on R'^ and by MsiW^) the set of 
signed measures on M*^. We denote by 5(Ms(M'^)) the set of symmetric matrices whose 
components belong to Ms(M'^). 

2. Framework and main results 

We consider an integer N ^ 2 and c g] — 1,;^y^]. We introduce a probability space 
(rj, P) and denote expectation by E. On this space, we consider a symmetric random 
matrix- valued Gaussian process X'^(x) = {Xfj{x))i ^ ij ^ n indexed by x G M'^ where the 
covariance structure is given for e > by: 

E[Xl,ixf] = j\ln - + 1), E[Xl,ix)Xl^ix)] = -cj\ln - + 1), i / j 

and for |y — x| > e: 

E[Xl,ix)Xl,{y)] =K{x-y), 
E[Xl,{x)Xl^iy)] = - cK{x -y),i^ j 

for some constant L > and some kernel of positive type K{x) = 7^ ln_|_ + g{x) where 
g is some continuous bounded function. In the sequel, we set g{0) = m. 

We also set a1 = 7^(ln j + ~ ^{y~^)- We assume that the off diagonal terms 

are independent of the diagonal terms and mutually independent of variance = ^^"^^ 
and covariance for |y — 3;| > e: 

E[Xl^ix)Xl^iy)] = ^K{x -y),i< j. 

In fact, the above structure is the most general situation that ensures that, for a given 
X € M"^, the Gaussian random matrix X'^{x) is isotropic (see Lemma 4 in the Appendix). 

Remark. Notice that the diagonal terms are independent if and only if c = 0. In this 
case, the above structure coincides with the usual Gaussian Orthogonal Ensemble (GOE) 
[1, 20]. Notice also that the boundary case c = jyi^ corresponds to trace-free matrices. 
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Remark. The canonical example of such a kernel K is when it coincides with 7^ In 
for X small enough. In dimension 1 and 2 we can even choose K{x) = 7^1n_|_i^. In 
dimension greater than 3, we can use the constructions developed in [15, 23]: for examples 
of such kernels, see Appendix A.l. Another approach is to use the convolution techniques 
developed in [24]. This does not exactly fall into the framework set out above because the 
convoluted kernel depends on e at all scales, i.e. for |x — y| > e. Nevertheless, this has a 
non significant influence on the forthcoming computations so that we also claim that our 
results remain valid for such regularization procedures. 



We want to study the convergence of the following random variable which lives in S{M, 

M'{A) = - I e^'^^Ux, A C (1) 

Ce J A 

where is chosen such that E[M''(A)] = where |^| is the Lebesgue measure of A. 

We will prove that the normalization constant has the following explicit form: 



1 r(i/2) 
ivr(iv/2) 
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Theorem 1. Let < 7^ < d. Then there exists a random matrix measure M which lives 
in S{Ms{R''-)) and such that for all bounded AcR'^: 

E[tr(M'(A) - M(A)f] 0. 
We also have the following asymptotic structure: 



with Vn = ^ ^ ^J^- P'^r^h ermore, we get the following equivalent for k'^ 2: 

\iiE[trM{B{id,i)f] 



In i ^->o 

where Q{k) = dk — 7^ ^(^""^) ^ 



m (3) 



Notice that it would be interesting to prove that this matrix-valued Gaussian multi- 
plicative chaos admits a phase transition as in the scalar case, which is likely to occur at 
72 = 2d. 

Remark. Application in turbulence. In the paper [6], the authors consider the following 
boundary case as a building block of their random velocity fields: 

2 8 2 ,T o 1 1 
7 = -X , N = 6, c = — = - 

'3 iV- 1 2 

where is found to fit experimental data for ~ 0.025 [6, 7]. Here, the zero trace 
property is reminiscent of the incompressibility condition imposed on velocity fields. 
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Conjecture 2. The power law spectrum of M is given by the following expression: for 
allq<^,yK cB{0,L),yie (0,1], 

E[trM{B{0,l)y] ~ Cgf^'^\-lni)^, 
where Cq > is a constant and the structure exponent is given by 




If this conjecture is true, this would show that non-commutativity yields an extra log factor 
in the power-law spectrum of M . 

Remark. Notice that one can define a notion of " metric" (actually a measure) through the 
quantity 

A G B{^'^) ^ tvM{A). 

Therefore we can define the notion of Hausdorff dimension asociated to this "metric" (see 
[11, 22]). It would be interesting to prove a corresponding KPZ formula and relate it with 
a KPZ framework . 



3. Proofs of the A^-dimensional case 



Let us first mention that several results abouts isotropic matrices and related computations 
are gathered in the appendix and will be used throughout this section. 



3.1 Joint law of the eigenvalues of Gaussian isotropic matrices 

We consider a symmetrical random matrix X = {Xij)i <^ ij <^ n niade up of centered 
Gaussian variables with the following covariance structure: the off-diagonal terms {Xij)i^j 
are i.i.d. with variance cj^. The diagonal term {Xi^i, • • • ,Xj\[^n) is independent from the 
off-diagonal and it has the following covariance structure: 

Kn = {E[Xi^iXjj])i ^ ij ^ ^ = (1 + c)ajl]y - cajPN 

where In is the identity matrix, P/v = (l)i,j and c e] — 1, j^l- By noting that = NPn, 
we get the following inverse for K if c ^ Jf^' 



""^^ ~ + c)^^ + aj{l + c) (1 + c{l - N)f^ 



The density of the random matrix, with respect to the Lebesgue measure {dxij)i <^ j, is 
therefore given by: 



1 1 V^JV ^.2 c 1 /T^N \2_ 1 „2 

f((x- ■) = —e 2<T^(i+c) 5^|(TT^IT+^(T^^'>'^«=i^'''-' -^^i<j^i 

where 



Zj, = (27r)^(^+i)/V^a^(^-i)/2(l + c)(^-i)/2yr3|iV^ 
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is a normalization constant. 

Therefore if we have the following condition: 

al{l + c) = 2a\ (4) 

as we have required in section 2, we can rewrite the above density in the following matrix 
form: 

1 7^ , ,\ (trX)^ trX2 

/((^M)i^.) = ^e (5) 

with Zn = 2^/27r^(^+i)/V^(^+i)/2(i + c)(^-i)(^+2)/Vl + c(l-A^). This shows that 
the matrix is isotropic, namely that for any real orthogonal matrix O, the matrices X and 
OX^O have the same probability law. Therefore by applying [1, Proposition 4.1.1, page 
188], we get the density of the unordered eigenvalues: 

f{{h)i^i^N) = ''^^(1+-) ni<j|Aj-Ai|, (6) 

^here a = ,,.(1^,) and = 2^(^-i)/^^^i^Jg)gf (notations of [1]). We 

remind that p{Un{^)) = 2^/2(2^)^(^+1)/^ HaLi (^ee [1, page 198]) and thus: 

Z^ = N\{2.ri\{\ + c)(^-^)(^+2)/Vl + c(l-iV). (7) 



k=l 



The isotropic condition (Eq. 4) ensures also that the collection of eigenvectors {vi)i ^ i ^ n 
is independent of the eigenvalues {Xi)i ^i^N-, and they are distributed uniformly on the 
unit sphere according to the Haar measure [1, Corollary 2.5.4, page 53]. 

3.2 Computations of the renormalization 

We consider here isotropic symmetric matrices X'^{x) = (X|j(x))i ^ jj- ^ at as defined in 
section 2 and compute the renormalization of order 1, i.e. the constant such that: 

The isotropic nature of the matrices ensures the proportionality of the former expectation 
to the identity matrix /^r . We want more precisely an equivalent of Ce as e — )• 0. We have: 



'N JrN 



where = 2a'^{i+c) (i+c(i-Af)) ^^^'^ normalization constant Zat given by Eq. 7 with 
We set Ui = — and therefore we get: 



(J. 

Ce = 



0"e 

N{N+l)/2 



Zn 
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where a = 2(i+c) (i+c(i-Af)) • introduce 



^{ux,--- ,un) = (T,ui - a{J2uif - + c) ^ 



i=l 

The function is maximal for ui = «S'e(l + 2a(l + c)(A^ — 1)), i ^ 2 : Ui = — 2aS'e(l + c) with 
•S'e = ^% AT - We thus set ui = ui + S'e(l + 2a(l + c)(A^-l)), i ^ 2 : = t;i-2aS',(l + c) 
to get: 

N(N+l)/2 fi „ „ 

X ^ j<jbi - Vi\dvi ■ ■ ■ dVM- 
Therefore, we get the following equivalent by using the Laplace method: 

JV(iV+l)/2. .jv-l iV-l 4 r ,2 1 v^iV 2 



I " — I J 

e-!>0 Z^v JK^ ^ ^ 

By using equation (25) in the appendix, this leads finally to the following equivalent as 
e ^ 0: 



1 r(i/2) 



>o ivr(iv/2) 



(l + c)(^-i)/Vf-ie^. (8) 



3.3 Computation of the moment of order 2 

In order to study the convergence, for e — )• 0, of the Gaussian chaos M'^{A) (Eq. 1), we 
need to consider first the second-order moment: 

E{M'{Af) = \ [ ^(e^^(^)e^'(^))dxdy, 
ci JaxA 

that involves the following quantity: 

In- (9) 

We will show that E(M''{A)'^) converges to a limit as e — ?■ 0. From this convergence, 
one can easily deduce that the sequence {M'^{A))^yQ is a Cauchy sequence. Again, the 
proportionality to the identity matrix in (9) comes from the isotropic character of matrices 
and we will see moreover that, because the so-defined field of matrices is homogeneous, the 
former quantity will depend only on |x — y| . The purpose of this section is to compute this 
quantity. We will restrict to the case |?/ — 3;| > e as the case |y — x| ^ e, once integrated, 
leads to vanishing terms in the limit e — )■ 0. It requires first the derivation of the joint 
density of the two matrices X^(x) and We will see indeed that the quantity will 

depend only on |a; — y|. We will also notice that, contrary to the one-point density (Eq. 5) 
from which it can be shown that eigenvectors and eigenvalues are independent, eigenvalues 
at point X are not only correlated to eigenvalues at point y, but also with eigenvectors at 



i^(e^^(-)e^^to)) = tr(e 
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point y. This intricate correlation structure is reminiscent of the non-commutative nature 
of this field of matrices and is encoded in the so-called Harish-Chandra-Itzykson-Zuber 
integral over the orthogonal group, or angular-matrix integral, and its related moments. 
This is an active field of research in random matrix theory and up to now, no explicit 
formula are known in dimension N ^ 3 (see for instance [5, 4, 8] and references therein). 
Nonetheless, we will succeed to get an explicit result in the asymptotic limit e — )• 0. 



3.3.1 Joint density of two isotropic matrices 

We consider here two isotropic symmetric matrices X^(x) = {X^j(x))i ^ ij ^ n and X^{y) = 
{X^ ,j{y))l ^ ij ^ AT as defined in section 2. We recall that matrix components are logarith- 
mically correlated over space. We note Xij = X^j{x) and yij = X^ ,j{y), and in matrix 
form X = X'{x) and Y = X'{y). 

Let us first consider the diagonal terms 

- ,XN,N,yi,i, - ■ ■ ,yN,N)- 

The covariance structure K2N of these elements is given by: 



K2N 



where Aj^ = (1 + c)In — cPn and we recall that = 7 (In 7 + 1) and = 7 In \ x-y\ ■ 

We know that the inverse of K2N is given by: 

7^-1 _ e N \y-x\ N 

where = jj^In + 2aPN with a = 2{i+c) (i+c(i-Af)) '^^ich leads to the following 
density: 

fiiXi,i)l ^i^N] {yj,j)l ^ i ^ Af) = 

_ <^?/(l + c) Ei ^|,i + 2a<T?{Ei ^i,«)^+<^?/(l + =) Ei vli + ^'^'^ciT.i !/i,i)^-2a2^_^|/{l + c) ^i,»!/i,»-4<'fy_^|«{E, ^»,i)(Ei y^,i) 

2(4-'^? 7) 

Cive if-^l 

= i2.r^LiK,^y det(i^27v) = {af - + c)2(^-i)(l + c(l - iV))2 

and therefore cat = (2^)iV(^4_^4 _ piv/2(\+e)(iv-i)(i+e(i-jV)) ' ^ similar procedure can be 
performed for the remaining A^(A^ — 1) off-diagonal terms of the two matrices. The density 
of the couple {X = X'^^x), Y = X^{y)) is thus given by, in matrix form: 



f{X,Y) = 

(10) 

where cat = cat ^jv(iv-i)/2Q , eW(jv-i)"/2(^4_^4 Njv(iv-i)/4 ■ We can see in the expression of 

the joint density of the two matrices X and Y (Eq. 10) two different contributions. The 
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first one, where is entering terms of the form tiX'^ + try^ and (trX)^ + (trY)'^, relates 
the density of two symmetric isotropic matrices as if they were independent. The second 
contribution relates an interaction term coming from the logarithmic correlation of the 



components. Indeed, the former vanishes if the matrices are independent, i.e. a\ 



y-x\ 



0. 



At this stage, it is convenient to introduce two i.i.d. random matrices M = (Mij) and 
M' = {M[ ■). These random matrices are taken to be living in the Gaussian Orthogonal 
Ensemble (GOE), namely they are symmetrical and isotropic with independent compo- 
nents with the following distribution: the components (Mjj)j ^ j are independent centered 
Gaussian variables with the following variances: 



E[Mi 



2a2 



-, i < j; E[MfA 



With this, we get the following expression for E[F{X{x), X{y))], where F is any functional 
of the two matrices X(x) and X{y): 



E[F{X{x),X{y))] 

OLO 

F{M,M')e *^ 



(l + c)(<7|-<Tf J 



tvMM'- 



trMtrM' 



where 



E 



((trM)2 + (trM')^)+ 



(l + c)(£r| 



\v- 



-txMM' 



i) 



-trMtxM' 



•:\ 



fill 



By using classical theorems about isotropic matrices, we know that M = OD{XfO, M' 
O' D{\'YO' where O (resp. O') is uniformly distributed on the orthogonal group of 
and is independent of the diagonal matrix D{X) (resp. D{X')) the diagonal entries of which 
are the eigenvalues of M (resp. M'). 



3.3.2 Joint density of eigenvalues of tw^o correlated isotropic matrices 

We are interested here in computing the renormalization constant Z (Eq. 11). To do so, 
we diagonalize the matrices M and M' , and perform an integration over the remaining 
degrees of freedom left by the eigenvectors (see [4] for instance) . We define the eigenvalues 
of M as A = (Ai, . . . , Aiv) £ IR'^ and we note the Vandermonde determinant as A(A) = 



^ i<j s£ N 



|Ai — Aj|. We get: 



Z = ^ [ |A(A)||A(A')|e 



2(1+c)(<t| 



2(l + c)(a« 



\y-^\' 



X e 



J{D{X),D{X'))d\d\', 



where i?^ is a renormalization constant such that 



R 



|A(A)||A(A')|e 



2(l+c)(<Tf 



Y-iV ,2 



2(1+c)(<t|-<t 



\y-^\' 



V^iV X /2 



dXdX' 



1, 
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and J is the following Harish-Chandra-Itzykson-Zuber integral [5, 4, 8], also called matrix 
angular integral {dO stands for the Haar measure on OAr(M)) 



r-JF4 rtr D{\)OD(X')0-^ 

J{D{X),Di\'))= I e'^+^""^-^i^-^l' dO, 



obtained while integrating over the eigenvectors that enter in the term tiMM' of Eq. 
11. We make the following change of variables Ui = , "\ Xi, u[ = i K (set 

7, = and get to: 

X e-"«^-^ -^^^^^-^ <f )^^^l. -)(E£. Diu'))dudu', 
where we have set: 

JiDiu),D{u'))= [ giTS^E5=i«.«;io,,,f 
iOjv(M) 

Therefore, since J{D{u),D{u')) converges pointwise towards as e — >• 0, we can use the 
Lebesgue theorem to get the following equivalent as e — )■ 0: 



Z ~ / |A(u)||A(n')|e"W^)^-i"-2(T+^2^.=i«i 

^ e-<ii:t.u.?H^t,<?)dudu\ 
that is straightforward to compute (see the appendix). 

3.3.3 Tvi^o-points correlation structure of the matrix chaos 

We want to get an equivalent as e — t- of the quantity given in Eq. 9. To do so, we 
consider the following quantity: 

r^^l ((trA-/)2 + (trM')2)+ trMM'+ '^r^V^ tvMtvM' 

Z = E[ti{e e )e If-^l Is'-^l l^'-^l J. 

In the same spirit as formerly, we diagonalize the matrices M and M' and perform the 
integration over the eigenvectors. We get: 

Z=\ |A(A)||A(A')|e ' ^oT^k^F^ ^'^^ » 

xe^^'i!/-!' I{D{X),D{X'))dXd\', 
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where / is the fohowing moment of the angular integrah 



■2 



r ^-4r tr D(X)OD{\')0-^ 

I{D{X),D{X'))= tr(e^WOe^(^)0-i)e<^+^'"'^-''l-^!' dO. 

JonW 

We make the fohowing change of variables Uj = . A^, n' = . A ■ (set 



_ V " \y-^\ - 



4_ 4 " * / 4_ , 



N NiN+1) 

Z = V / |A(n)||A(n')|e"W^)^fc=i"fc-2(T+^^fe=i"* 



-"((Er=i"^)'+(E^=i<)')+ 
where we have set: 



2 



r 1 \y—x\ I 1^ |2 

JOjv(K) 

known as the Morozov moment [4]. We make the following change of variables in the 
above integral: Uj = + 7^, Uk = Vk — cje for / i and u'j = v'j + 7e, u'j^ = v'j^ — C7e for 
k ^ j. We obtain the following equivalent: 

^ ~ E / |A.(^)I|A,K)I 

where |Ai(v)| = 1^' ~ ^''1 

Jonim 



af .(ilL-2c) f 1^ |2 fTp i(l+c)|Oi iP ,^ 



which is independent of Therefore, we get: 

7V(7V+1) 2 ^2{iV-l) 2(iV-l) 

Z ~ iV^^^ ^il+i^^ V / |Ai(.)||AiK)| 

Cat Jm^^xirJv 

g- 2(1^^ ELl -i-5(TT^ ELl <'-"((ELi -fe)' + (ELi .D^)+2aaf^_,(l+c(l-JV))2^^^^,_ 

In conclusion, we get: 

r(A'/2) JojvCR) 
Including furthermore the normalization constant Ce (Eq. 8), we get: 

Z/{Zcl) ~ iV2e— Ui /■ |Oii|2e'^l.-|(i+^)|Oi-^l'dO. 

iO]v(M) 
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3.3.4 Computation of the moment of order 2 

From the above subsections, we deduce that: 

J Ax A yOiv(M) 

We remind that the law of |Oi^ip is the one of the square of one component of a vector 
uniformely distributed on the unit sphere, and has thus a density given by (see Lemma 3) 

f(v) = v-'/Hl - v)^''-'''>/^. 

r(i/2)r((iv-i)/2) ^ ' 

We get the foUowing equivalent as |y — x| — )■ 0: 



jo^m \y--\^^ r(i/2) (i + c)(^-i)/2a^ri ' 

which entails (2). 

3.4 Computation of the moment of order k 

We are interested here in studying the convergence, when e — )• 0, of the Gaussian chaos 
M''{A) (Eq. 1) for higher order moments such as, A; E N, 



^X^Xi) j . . . ^^^^ 

k 



that involves the following quantity: 



4 n 



,1 < i s£ 



tr JJ e 



1 



In. (12) 



To generalize former calculations in the case k = 2, we will first derive the joint density 
of /c-matrices (X^(xj))i ^ j ^ fc. A generalized version to A;-points of the Harish-Chandra- 
Itzykson-Zuber integral enters the expression of the density. An exact evaluation of these 
integrals remains an open issue. As far as we know, only their behavior in the asymptotic 
limit of large matrices (A^ — )• +oo) has been considered in the literature [8]. Nonetheless, 
a logarithmic equivalent of the quantity of interest (Eq. 12) can be obtained and allows 
to show the multifractal behavior of the multiplicative chaos (i.e. C{k) is a non linear 
function of the order k, see theorem 1). 

3.4.1 Joint density of k isotropic Gaussian matrices 

We consider here k isotropic Gaussian matrices (X'^(xi))^ <i^ k- '^^^ ensemble made of 
the kN diagonal terms, i.e. 

{Xli{xi),--- ,X^^^(xi),--- ,Xli{xk),--- ,Xlf,N{xk)), 
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has covariance structure i^fciv^ 



KkN = 













a^AN 




'^|X2-X,1^^ 











N 



'^\Xk-Xk-l\'^^ 



N 



J 



where again, An = [1 + c)In — cPn- We know that the inverse of K^^ is approximately 
given by (e — >■ 0): 



<^\x2-xi\^N 



-1 



'^\X1-X2\^N 

Oe^N 



— CTi 



v 



|a;fc_i-a;i| N 
^\xk-xi\^N 



A J -a\ 



A' 



-1 

\xk-l-X2\^N 

^|xfc-a;2riV 



^2 



n-24-1 
Nfc-a:fc-i| ^ 



\X2-Xk\ N 



-rP- 4-1 



where ^jv^ = j^^^N + 2aPAr, with a = 2(1%) (i+c(i-jv)) • '^^^ density of diagonal 
components, considering the AT-dimensional vector X^^ = is 
thus given by: 



where cjv 



)*^'''(XTT^^JV+2aPAr)xO) 

For the off diagonal terms, the situation is simpler. If 



(27r)'=^/Vdct(i^fciv) 

i < j, the covariance matrix of the vector {Xf ,j{xi), ■ ■ ■ ,Xf j{xk)), which is independent 
on all the remaining diagonal and off-diagonal components, is: 



\X2—Xl\ 



\xi-X2\ 



'''|xfc_i-a;i| '^\xk-i-X2\ 

\ '^\Xk-Xl\ '^\Xk-X2\ 



cr: 



a 



\Xk-Xk-l\ 



\xk-l-Xk\ 



whose inverse is approximately given by (e — >■ 0) 



—a 



\X2-Xl\ 



— cr I 

\Xk-l-Xl\ 



\x1-x2l 
-I 



—a 



\Xk-Xk-\\ 



\X\-Xk\ 



\Xk-\-Xk\ 



This leads to the following density, using the notations xf^ = Xf^(xr): 



/( 
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Therefore, we get the following density for the k matrices (we omit superscript e for the 
sake of clarity) : 

which we rewrite under matrix notation: 

f{X{Xi),- ■ ■ ,X{Xkj) = CNC 2(1+=)-. 

We introduce k i.i.d. random matrices M^'^ = (M^'j) pertaining to the GOE ensemble. 
These random matrices are symmetrical and isotropic with independent components with 
the following distribution: the components {M^ 'j )i ^ j are independent centered Gaussian 
variables with the following variances: 

= ^ < r, EmSf] = (1 + c)ai 

With this, we get the following expression for the expectation of any functional F{X^{xi), • • • , X^{xk 
of the A; matrices X^{xi), - ■ ■ , X'^{xk)- 

E[F{X'(x^),--- ,X'{xk))] 

" ' ' Z ' 

where: 

By using classical theorems about isotropic matrices, we know that, for each r, M^^) = 
O^^"^ D{X^^'^yO^^^ where O*-''-' is uniformly distributed on the orthogonal group of and 
is independent of the diagonal matrix D(X^^^) the diagonal entries of which are the eigen- 
values of mW. 

3.4.2 Joint density of eigenvalues of A; isotropic Gaussian matrices and com- 
putation of the renormalization 

We start by computing Z. For A = (Ai, . . . , A;v) € W^, we note the Vandermonde deter- 
minant as A(A) = Hi ^ i<j ^N\^i~ ■^jl- We get: 
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where i?^ is a renormalization constant such that 



ne / -r=ii — V- /i- -^J 



and J is the following angular integral: {dO^^^ stands for the Haar measure on On{M.)) 
J(D(AW),--- ,I5(aW)) = 

Jon{K)'^ 

(r) A^''-* 

We make the following change of variables ul = 

N(N+l)k/2 r 1 ^fc s^N ,M.2 „Y^fc (-r^N (r).2 



Z = ^^-^ / n,^=i|A(uW)|e-TO 



where we have set: 

J(L>(u«),--- ,Z?(u(*^))) = 

Therefore, since J(£)(u'^^)), • • • ,D{u^^^)) converges pointwise towards 1 as e — >■ 0, we can 
use the Lebesgue theorem to get the following equivalent as e ^ 0: 

NiN+l)k/2 1 , W., 

3.4.3 fc-points correlation structure of the multiplicative chaos 

For i ^ j, we want to get an equivalent as e — >^ of the following quantity: 

We get: 

xe^^'^<'"i--^<i^^^=i^' ^^^^=^^» ^/(Z)(A(i)),--- ,L>(A('=)))dA(i)---dA('=), 
where / is the following angular integral: 

/p(AW),--- ,D(AW)) = / (n.^^iOWe^^^^'"^)*©^),,,- 

JO]v(R)'= 
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(r) A^''' 

We make the following change of variables = 
'o=»;'fc=J 



where we have set: 



Wc make the following change of variables in the above integral for 1 ^ r ^ k: 

(r) (r) (r) 

VjJ + CTg, Uk = — cae k / jr- We obtain the following equivalent: 



Z ~ 

e->0 



N JV(JV+l)fc/2 fe(,2/2.. x(jV-l)fe^(iV-l)fe 



^ g- 2(tW 2:^=1 Ef=i{-r')^-aE^=i(Ef=i .r)^+2a(l+c(l-iV))2 ^^^^ af^^^^. 



AT 

X 



where wc have set: 



pi,---,jk _ / ('jjA; 

lo,h,— ,lk-l,lk J (n\k lr-l,jr Ir ,. 



X e 

sir Ji 



Oiv(R)* 

1 ^2 v-Af ajr,jl spN ^(r) ^(r) ^(l) ^(i) 

(1+c) Z^r-<! °^|xT--a:(| ^iTii,m2=l Pm,i,m2 Z^ni,n2 = l '-'ni .mj l-'n2,mi l-^nj ,m2 <-^n2 .mj (JQll) . . . dQW 



with = (-C + (1 + c)l^i=jj(-c + (1 + c)lm2=h)- This leads to 



jjl-!'" fjk 



[ (ntio;:!,,o;:)j 

^g(li^E.<i<^f._.,|(c^*K0«*0W0«*0W)+(l+c)2(*0W0W)2^_.p 
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In conclusion, we get the following equivalent: 
where we have the following expression for Fij{xi, ■ ■ ■ , x^): 

N N 



We get thus the following expression: 

Z/(Zc,^) ~ Cfc,^e"'^^<'"l---^IF,j(xi,--- ,Xfc), (13) 

where C^^tv is a constant which depends only on /c, (we can compute this constant but 
it is tedious and will not be necessary for the purpose of this paper). 

3.4.4 Computation of the moment of order k and of the structure functions 

From relation (13), we get the following expression: 

E{trM'{Af) ^ Ck,N e""^'-<'"l---H \^F,,,(xi,--- ,Xfc)dxi---dxfc (14) 
The main difficulty is to study the functions Fjj. If one takes the trace, we get: 

N 



^ ^ ) X/j) 

E / (*o«o(^)),,,, . . . (*o('=-i)oW),,_,,,(*o(^)o(i)), 



1=1 

N 



X e(^+^)^-'<-^.|('^^^^«'')->^.dO(^) • • • dO^'l (15) 

In particular, for k = 2, we recover that: 

N „ 

Vf,,(xi,X2) = iVM (*0«0(2))2^e(^+^W..-..|(*0*^'0''')?ad0Wd02. 



i=l 
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In order to prove (3), we want to study the asymptotic of Yl^=i ^iii^^ii ' ' ' A^k) as 



with (xi, • • • , Xk) fixed. Since Y.vA^O'^^^O^^'^)'] i < ^^^> it is obvious that: 



e^o In i ^ ^ ^ ' 2 ^ ^ 

To get an inverse inequahty in (16), we will study the asymptotic of each term in the 
sum (15). Here we suppose that L = 1 and m = to simplify the presentation. We fix 
(ill ■ ■ ■ 5 Jfc) and e, 6 small such that e < 5. We have: 

(*o«o(^)),,,, • • • (*o('^-i)o('^)),,_,,,(*oWo(i)),,,, 

= A, + A,^s + As 

where 

A= ■■■ , A.s -- , 

and is the l^r<«(*0^''^0*"''')ir,j; ^ M^ili _ § part of the integral. On the event 
Er<iCO^''^0'-^^)lj^ ^ - e, each I^O^'^^OC))^, jj is greater or equal to VT^. In 

particular, we have that | (*0('^)0('-+i))j, j^+J ^ \/l — e for all r ^ k — 1. Notice that 
(tO(fc)0(i))j-^j^ = (*0(i)0(2))j-^_j2...(*0(^-i)0W)j^_^j^ +0(e). Therefore, we can con- 
clude that (*0(i)0(2))^.^^^_..(to(fc-i)o(fc))^.^_^^^.^ j^j^^i (*0(i)0(2))^-^_^-2 have the same sign. 
Thus, we get: 

A, 

(*o(i)o(^)),,,v(*o('^-^)o(^)),,_,,,(*o('=)o«),,,, 

■>'E.<,(*oMo«)?,,,,, > 

The only thing to check is that X]r<«(*^*'''^^^'^)ir ji ^ ^^^2~^^ ~ ^ ^^^^ ^ positive probability 
but this can be seen easily by setting one chosen element of each 0^"^^ say , very 
close to one. Now, one can choose 6 larger than e such that \A^^s\ ^ 4^- Finally, for these 
choices of e,5, we have: 
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We thus get the following: 

r HAe + A,,s + As) ^ k{k-l) 
hm ^ (1 + c)7 ( e) 

Since this is valid for all e, by combining with (16), we get that: 

™ r~i = + — 5 — ^^^^ 

^-s-o In J 2 

The desired result (3) is then a consequence of (17) and of the limit (14). 

A. Appendix 
A.l Discussion about the construction of kernels 

In this subsection, we discuss in further details the construction of the kernel K as sum- 
marized in remark 2. In dimension 1 and 2, it is plain to see that 

ln+— = / {t-\x\)+ULidt) (18) 
1^1 Jo 

where the measure vl is given by (5„ stands for the Dirac mass at u): 

Hence, for every ^ > 0, we have: 

L 1 

ln+ — = -ln+ — - = / (t - |x|^)+z^LM(dt). 
\x\ /i \x\^^ Jo 

We are therefore led to consider fi > such that the function x i— t- (1 — |x|'')_|_ is positive 
definite, the so-called Kuttner-Golubov problem (see [14]). 

For d = 1, it is straightforward to check that (1 — |x|)+ is positive definite. We can 
thus consider a Gaussian process X"^ with covariance kernel given by 

fL 

.2 



K,{x) =Y I {t- \x\)+iyL{dt). 

Notice that 



and 

Ve 



Vx / 0, 7^ ln+ = lim K,{x) (19) 

\x\ e-S-O 

<\x\^L, K,{x)=j^ [ it-\x\)+i^Lidt)=jHn+^. (20) 

J\x\ \X\ 

In dimension 2, we can use the same strategy since Pasenchenko [21] proved that the 
mapping x i— )■ (1 — |x|^/^)+ is positive definite over R-^. We can thus consider a Gaussian 
process with covariance kernel given by 

/.Li/2 

K,ix) = 2j^ (t-|x|V2)^z,^,/,(dt), 

Jel/2 
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sharing the same properties (22) and (23). 

In dimension 3, it is not known whether the mapping x i— t- ln+ admits an integral 
representation of the type explained above. Nevertheless it is positive definite so that we 
can use the convolution techniques developed in [24]. In dimension 4, it is not positive 
definite [24] so that another construction is required. We explain the methods in [23]. We 
set the dimension d to be larger than d ^ 3. Let us denote by S the sphere of M.'^ and 
a the surface measure on the sphere such that cr(5) = 1. Remind that this measure is 
invariant under rotations. We define the function 

Va; G M.'^, F{x) = 7^ / ln+ --^^a{ds). (21) 

Js iK^^Vl 

It is plain to see that F is an isotropic function. Let us compute it over a neighborhood 
of 0: for |x| ^ L, we can write x = \x\ex with Cx G S. Then we have 

f L L f 1 
F(x) = 7^ / In—r-- — a(ds) = In — + / In-- —a(ds). 

Js m\{ex,s)\ \x\ Js \{ex,s)\ 

Notice that the integral In ^^^^-^^(ds) is finite (use Lemma 3 below for instance) and 
does not depend on x by invariance under rotations of the measure a. By using the 
decomposition (18), we can thus consider a Gaussian process X*^ with covariance kernel 
given by 

„2 



K,{x)=Y I I {t-\{x,s)\)+UL{dt)aids), 

sharing the properties: 



'S 



and 



Vx / 0, lim K,{x) = F{x) (22) 

Ve<|x|^L, KJx) = F(x) = X^ln^ + C (23) 

\x\ 



for some constant C. 



A. 2 Auxiliary results 

We give a proof of the following standard result: 

Lemma 3. If {Zi)i <g i ^ tv tJ'^e i.i.d. standard Gaussian random variables then the vector 

y= I ^ (Zi,...,Zjv) 

is distributed as the Haar measure on the sphere of . In particular, the density of the 
first entry of a random vector uniformly distributed on the sphere is given by: 

r(-^) _i , jv-a 

r(i)r(^) ^ 5(1-2/) 2 i[o,i](2/)dy. 
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Proof. By using the invariance under rotations of the law of the Gaussian vector (^j)i ^ i at, 
the law of V is invariant under rotations and is supported by the sphere. By uniqueness of 
the Haar measure, V is distributed as the Haar measure. We have to compute the density 
of Ci = „a; ^ ^1 ■ Notice that 

. _ y 
Y + Z 

where Y, Z are independent random variables with respective laws chi-squared distribu- 
tions of parameters 1 et — 1. Therefore 

E[f{Ci)]= I I /(^^)^^x-5e-f^^3^i— — yVe-f dxdy 
Jr+Jr+ ^x + yJ 22r{^) 2— r(^) 

If y N-2 



2Tr(i)r(^) Jo " ^ 'v^(i-n)f 

J{U)U 2 (1 — ti) 2 du. 



f{u) 3" / e 2(i-«)y 2 dydu 



r(i)r( 







□ 



Next we characterize all the symmetric Gaussian random matrices: 



Lemma 4. Let X he a centered Gaussian symmetric random matrix of size N x N. Then 
the diagonal terms {Xu, . . . ,Xnn) have a covariance matrix of the form (T^(H-c)lAr— ccr^P 
for some ^ and c g] — 1, jfrj]! where P is the N x N matrix whose all entries are 1. 
The off- diagonal terms are i.i.d with variance cr^^^^ and are independent of the diagonal 
terms. 

Proof. If X admits a density with respect to the Lebesgue measure dM over the set of 
symmetric matrices (see [1, chapter 4]), then the density of M is given by 

e-^(^-^) dM, 

where / is a homogeneous polynomial of degree 2. By isotropy, / must be a symmetric 
function of the eigenvalues of M. Therefore it takes on the form 

/(M) = atr(M2) + /3tr(M)2 

for some a, /3 G M. In this case, the result follows easily. 

If M does not admit a density with respect to the Lebesgue measure over the set of 
symmetric matrices, we can add an independent "small GOE", i.e. we consider M + eM' 
where M' is a matrix of the GOE ensemble with a normalized variance independent of 
M. The matrix M + eM' admits a density so that we can apply the above result. Then 
we pass to the limit as e — 0. □ 

A. 3 Some integral formulae 

Let a, c > 0. We want to compute the integral: 

^-»iI:tl^^f-2(T+^^^=l^^Ui<j\Xj - Xi\dX. 
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We write the integrand under the form (6): 

where (T^(1 + c) = (1 + c) and a = 2a^(i+c) {i+c{i-N)) • -'■'^ ^^st case, we have c = 
i+2a(°S(i^-i) ^'^^ 1 + c(l - iV) = i+2«(i+c)(jv-i) - We deduce: 

jR'v ^^^^ i Y^l + 2a(l + c)iV 

(24) 

We also want to compute the integrah 
We have: 

/ e~°^^'=i^^^'"W^^»=i^»' n |Aj - Ai|dAi---(iAiv 

= / ( /"e-2"Ai(Ef.2A0-(a+5^)A?^^^^g-a(Ef=2A,)^-2(I^Ef.2A2 | A^- - A^ | dAa • • • dA^ 

= V2^j / + 1/ e^"^^^^(^-^-)^-"(^-^')^-^^-^? n |A,-A.|dA2---dA.. 

Y 2a(l + c) + 1 iiRiv-i 2^i<i' J *l ^ 

= V2^J — -^-^^ / e" (^-2 A,)2-5(TVjy E^a A? n |A,--AJ(iA2---(iA^ 

Y 2a(l + c) + 1 V-i 2^i<i' J *l 2 

Y 2a(l + c) + 1 JiRiv-i 2 ^ i<j' 



for cj^(l + c) = 1 + c and c = 2a{i+})(T~i)+i (o^ equivalently, l + c(2-A^) = i+2j"(i°iyj^Li) 
and 1 + c = 1 4.2^^^! l-ex'jv - 1 ) ) • '^^is leads to the following: 

7{Eil2 Ai)^-:r:57TT~V E,l2 Af 



e d 



5^j^a+.(2-iv))Vz.,.2 i^j(T:^z.,=2 . n JA, - Ai|dA2 • • • dA^ = Z^v-i 



2 s: i<j 



where = (iV-l)!(2^)(^-i)/2(nf=Y n^)af ^-^)/^l+c-)(^-2)(^+i)/Vl + ^"(2 - iV). 

In conclusion, we get: 

/" e-°(E£iA,)^-2(iVj)E.=i^? n |A,- - Ai|dAi---dA7v 

Jrat 2 s: j<j 

^ ^ r(l/2)^^l + 2a(l + c)iV ^ ^ 
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